function test13(  )
%TEST4 Summary of this function goes here
%   Detailed explanation goes here
    
    
    W = imread('../TestImg/stormcolor.jpg');
    
    W = double(W(:,:,:));

    S = fspecial('gaussian',7,7);
    
    %H2 = imfilter(W,S);
    N = [size(W,1)+size(S,1)-1, size(W,2)+size(S,2)-1, size(W,3)];
    H2 = zeros(N);
    H2(:,:,1) = conv2(W(:,:,1), S);
    H2(:,:,2) = conv2(W(:,:,2), S);
    H2(:,:,3) = conv2(W(:,:,3), S);
    
    r = 0.9995+(rand(size(H2)).*(1.0005-0.995));
   
    H2 = H2.*r;
    
    t = cputime;
% %     WRec = MyDeconv2DFFT(H2,S);%.*sum(S(:));
%     WRec = zeros(size(W));
%     WRec(:,:,1) = MyDeconv2DFast(H2(:,:,1),S,size(W,1),size(W,2),1000,0);
%     WRec(:,:,2) = MyDeconv2DFast(H2(:,:,2),S,size(W,1),size(W,2),1000,0);
%     WRec(:,:,3) = MyDeconv2DFast(H2(:,:,3),S,size(W,1),size(W,2),1000,0);

%    WRec = MyDeconv2DImage(H2, S, 'RichardLucy', 10000);
   
   WRec = MyDeconv2DImage(H2, S, 'RichardLucy', 60000);

    disp(cputime-t);
    
    
    figure(1); 
    subplot(2,2,1);
    imshow(ZeroToOne(W));
    disp('W');
    %disp(W);
    
    subplot(2,2,2);
    imshow(ZeroToOne(S));
    disp('S');
    %disp(S);
    
    subplot(2,2,3);
    imshow(ZeroToOne(H2));
    disp('H');
%     disp(H);
    
    subplot(2,2,4);
    imshow(ZeroToOne(WRec));
    disp('WRec');
%     disp(WRec);
    
end

